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METHOD  OF  DETERMINING  ERROR  IN  CALCULATED 
PARAMETERS  OF  COMBUSTION  PROCESS  RESULTING 
FROM  ERROR  IN  THERMODYNAMIC  PROPERTIES  OF 
INDIVIDUAL  SUBSTANCES 


V.  Ye.  Alernnsnv,  A.  K.  Dregalln, 

Z.  Kn.  Gruzdeva,  and  A.  S.  Lyashev 

Kazan ' 


Thermodynamic  calculation  is  broadly  used  to  determine  the 
parameters  of  the  combustion  process.  Error  in  any  parameter  of 
the  process  4> , determined  as  a result  of  thermodynamic  calculation, 
depends  on: 

a)  the  mathematical  error  of  the  calculation  method; 

b)  error  in  the  original  data  on  the  fuel  (enthalpy  and 
elementary  composition  of  the  fuel); 


c)  error  from  reference  data  on  properties  of  individual 

substances  - components  of  combustion  products  (enthalpy  1^, 

0 ^ 
entropy  3^ , and  other  properties  of  q-th  substance). 


Studied  in  the  work  is  the  effect  of  error  from  reference 
data  on  the  properties  of  individual  substances;  errors  of  type 
a and  b have  already  been  discussed  in  the  literature  [1], 


In  the  case  of  a constant  elementary  fuel  composition  (number 
of  atoms  of  the  i-th  chemical  element  in  a molecule  of  the  fuel 
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r,^= ; -..ost } ar.y  parameter  of  the  combustion  process  4>  may  ivpreeen; 

the  function  of  enthalpy  i and  entropy  s of  the  combustion 

products.  The  enthalpy  and  entropy  in  turn  are  determined  by 

the  composition  of  the  mixture  and  by  the  thermodynamic  properties 

of  the  individual  substances  and  S^.  Consequently,  error  A# 

which  develops  as  a result  of  errors  in  the  thermodynamic 

properties  of  the  individual  substances  can,  in  the  final  analysis, 

be  expressed  as  AI^  and  AS^. 

q q 

The  reference  book  [2]  gives  the  following  expressions  for 
enthalpy  1^,  the  change  in  enthalpy  (H^-Hq)o,  and  entropy  of 
the  q-th  substance: 

/•-A //%,.,  + (//«  -//§), -(//^ -//$),;  (1) 

<*f (—'-)  . (2) 

... 

where  Q is  the  statistical  sum  with  respect  to  the  states  of 
q 0 

the  molecules;  AH  ^93  " the  heat  formation. 

As  we  see,  errors  AI^  and  AS^  are  determined  by  the  errors 

q q 

of  four  quantities,  which  will  henceforth  bo  conditionally 
denoted  as: 


V-Atfo/n,.,;  (H^  -//?>,;  In 


Consequently,  the  error  of  parameter  $ for  small  error  o> 

I I 

is  the  linear  functional  for  Aw  , 

qk 

= V V(-4-)  ■Ad).  f c \ 

l <s*  k »*'  (5) 

where  quantities  Au>ql(  are  determined  by  the  reference  data  of  [2]; 

- the  parameters  which  determine  the  form  of  the  combustion 
process . 


To  calculate  A<p  In  the  general  case  we  must  solve  the  problem 
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M'  linear  programming,  since  rev  certain  component  s AH  f,0, 
expressed  through  the  enthalpy  e:*  key  iiubr.t  nncer. , and  imto;’;'. 
A(hJ-Hq)  and  are  te  a certain  extent  interrelate.!. 


Let  us  obtain  some  general  relationships  for  the  derivatives 
O^/aw^)^.  As  we  know,  the  thermodynamic  description  of  many 
processes  in  closed  systems  is  based  on  the  assumption  that  a 
certain  function  of  state  tp  (for  example,  the  enthalpy  of  the  fuel 
for  isobaric  combustion,  entropy  - for  isentropic  expansion)  does 
not  change  during  the  process.  For  such  cares  any  parameter  of 
the  combustion  products  $ may  be  represented  by  the  functions 


f~-  i (p,  T . •„),  \ 

v Hp<  *•«»)  I 


(6) 


After  we  have  equated  the  total  differential  functions  of 
a and  6 for  p and  4>«const  we  get 


f-*U  + (A\  (_£L\  , 

\ <>",»  )p.  ♦ l «*••«»  Ip-  T \6T  Ip.  *f  \ <*“**  Ip.  + 


(?) 


On  the  other  hand,  if  we  write  the  total  differential  of  the 


function  ^(p,  T,  u^)  for  p and  t|>*eonst , we  get 


/ <>r  \ _ VS  Ip.  t 

\ 

Ur/,.., 


(8) 


Thus,  for  partial  derivatives  of  the  tnermodynamic  parameters 


the  following  expression 

is  valid 

l 

/ * \ 

_ t * ) _ 

l_£L\  . , 

(tH-. 

(9) 

\ <S»  Ip-  ♦ 

V Ip.  r 

\ * Ip-  T 

(«L) 

Let  us  use  formula  (9)  for  the  case  of  isobaric  combustion, 
described  by  the  equations 


pt  = const, 
/,  = const 


(10) 
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where  the  subscript  "k"  refers  to  parameters  in  the  combustion 
chamber. 


As  a result  wo  get: 

for  the  temperature  of  t.lie  combustion  products 

\ '*"«»  Ip.  i tp  V '*"«»  Ip.  t 

for  the  molecular  weight  of  the  combustion  products 


(id 


(ID 


for  the  number  of  moles  <d‘  the  q-th  substance  in  the 
combustion  products 


/ \ _ /jln.j  _ \_  / On,  \ 

\ Ip.  i v *•«»/*  r cp  \ r*-,f  )p.  r \ if.  Ip.  -f 


(13) 


where  a Is  the  isobarlc  expansion  coefficient;  c - heat 
P p 


capacity  at  p»const. 


The  derivatives  (3i/3w  . ) and  (3u,/3u>  . ) are  written  as  the 

qk  k qk 


result  of  differentiating  known  expressions  []].  For  example, 
for  enthalpy  we  liave 


y«.< (”-"•)  ...lsl) 

I « \ . -r  i 

* '*■»»  Ip.  r r,M,  \ )p.  r’ 


(ID 


where  M is  the  number  of  moles  in  the  fuel,  giving  us  the 
equality  of  [lj 


P I*,. 

/ 1>  In  i< 

\fl  In  x 

Ip  \ ''III  t 1, 

( i‘0 


where  u is’  the  molecular  weight  of  the  fuel;  x - any  parameter. 


To  determine  partial  derivatives  (3  In  n /3w  . ) _ and 

q qk  p , l 


(3  In  H ,/3w  . ) 


qk 


> T we  differentiate  the  system  of  equations  for 

P , 1 


thermodynamic  equilibrium  [1]  with  respect  to  at  p and  T«const 
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ed  equation  syst<  u ! r calculating  tin* 


U7) 

(18) 


where  6,^  and  are  the  Kronecker  symbols;  1 - atomic  component; 
,!  - molecular  component;  quantities  Kj  .and  are  determined  by 
dl  !'!'<-rent  lat  InK,  the  logarithm  >>i*  the  equl  1 l l»rl  urn  constant,  and 
ir>  •-’von  In  Table  1. 


Table  1 
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m-t  i kH**,  I 

1 1 

"*tm-  "l 

r 

«|c 

j 

r-n 
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V 
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~ *r 
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1 
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i 

e 

V* 

The  values  of  quantity  Aw.  an  be  selected  on  the  basis  '»* 

*t'  1 

put  ilahed  data.  :-'or  the  cal  •«.  aions  below  the  data  of  [2]  were 
a?ej  a • onr  with  materials  fri  m t he  ■ r i edit  * *n  -f  this  handbook, 
which  Is  belnr  prepared. 

We  took  errors  MAH^f^^  ) and  directly  from  the 

tabulated  data  of  [2];  error  In  W/N  was  determined  from  the  ratio 


(19) 


where  ••  Is  the  reduced  thermodynamic  potenlal,  whose  error  Is 

£iven  in  [2]. 

To  calculate  errors  b(d  In  we  used  the  approximate 

estimate  recommended  In  f2] 
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AS^(2  3)d'l>*#. 

Consequently,  the  greatest  error  equals 


(20) 


I 

l 


A 


(21) 


The  presented  method  of  estimating  the  errors  in  parameters 
was  used  to  determine  errors  AT^  and  Ay^,  which  arise  in 
calculating  the  combustion  process  of  several  fuels.  Characteristic 
quantities  for  relative  errors  57^2  and  Sy^J  are  given  in  Table  2 
(pressure  p,  - MN/m2 , a - excess  oxidant  ratio). 

k OK 


Table  2 


Fuel 

•m 

P» 

*r„.  % 

W * 

Oi+Hi 

nr. 

13 

037 

0.33 

N,0< + (CH.)iNNHi 

0,8 

i 

ts 

l.ll 

ov 

o.+kerosene 

0.8 

IS 

0.14 

0.78 

.j^ 

■ii 

;| 


a 


When  ^*s  the  given  method  can  also  be  used  to  estimate  errors 
in  the  calculated  parameters  of  isentropie  expansion,  particularly 
to  determine  error  in  specific  pulse. 
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